import numpy as np


a = np.array([1,2,3,5,2,3,4,1,5,6,7,3,2,1,4,4,6,7,3,2,4,7,8,100])
avg_vaule = np.average(a)
std_value = np.std(a)

# P(|x-u|>=3y)<=0.003
def abnormalData(x,u,y):
    if abs(x-u)>=3*y:
        return x
    return -1
for i in a:
    if abnormalData(i,np.average(a),np.std(a)) != -1:
        print(i)